Probing non-Abelian statistics of Majorana fermions in ultracold atomic superfluid 
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We propose an experiment to directly probe the non-Abelian statistics of Majorana fermions by 
braiding them in an s-wave superfluid of ultracold atoms. We show different orders of braiding 
operations give orthogonal output states that can be distinguished through Raman spectroscopy. 
Realization of Majorana states in an s-wave superfluid requires strong spin-orbital coupling and a 
controllable Zeeman field in the perpendicular direction. We present a simple laser configuration to 
generate the artificial spin-orbital coupling and the required Zeeman field in the dark state subspace. 
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Quantum statistics is a fundamental concept in physics 
that distinguishes fermions from bosons. The discov- 
ery of the fractional quantum Hall (FQH) effects opens 
up the possibility to study exotic quasi-particles beyond 
bosons and fermions, which are called anyons [J 
There are two types of anyons, with Abelian or non- 
Abelian statistics. For Abelian anyons, the winding of 
one quasi-particle around another gives a phase factor to 
the wave function which is a fraction of 27r. For Non- 
Abelian anyons, a winding of one quasi-particle around 
another is associated with a unitary transformation in the 
subspace of degenerate ground states [3]. Such unitary 
transformations form a representation of the braid group 
where the elements in general are non-commutative. Re- 
cently, the systems with non-Abelian anyons have at- 
tracted considerable interest as candidates for realization 
of fault-tolerant quantum computation 0, [Hj . 

Probing quantum statistics is a challenge for both 
Abelian and non-Abelian anyons. Recent theoretical 
works suggest that interference experiments in a well de- 
signed point contact geometry may serve as a probe of 
Abelian Q or non-Abelian anyons 0] in the FQH system. 
Besides of the quasi-particles in the FQH system, a vortex 
excitation in + ipy superconductor traps a zero-energy 
bound state which represents a kind of self-conjugate par- 
ticles called Majorana fermions [s^]. Majorana fermions 
have the feature of non-Abelian statistics and have raised 
significant interest in condensed matter physics [1,|1-12|. 
Several inspiring proposals have been suggested very re- 
cently to realize Majorana fermions in different materials 
lol 11 1 . A scheme to observe Majorana fermions has also 
been proposed in the p-wave superfluid state of ultracold 
atoms [121] . 

In this paper, we propose an experiment to probe non- 
Abelian statistics of the Majorana fermions with an ul- 
tracold atomic gas in the s-wave superfluid state. In con- 
tract to the p-wave superfluid state, which still remains 
to be realized experimentally due to challenge associated 
with instability of the p-wave Feshbach resonance fl3| . 



the s-wave superfluid state is stable and has been ob- 
served in many experiments [l^. It has been noted re- 
cently that an s-wave superfluid with spin-orbital (SO) 
coupling emulates properties of the p-wave superfluid 



[lol 11, 15, 16 1 . Realization of Majorana fermions in 
an s-wave superfluid state requires SO coupling and a 
controllable Zeeman field on the effective spins ^11]. For 
ultracold atoms, an effective SO coupling can be gener- 
ated through gradient of spatially varying Berry phase in- 
duced from laser beams ll7Hl9| . We propose a scheme to 
realize both SO coupling and the required perpendicular 
Zeeman field on the effective spin, which a laser configu- 
ration compatible with the current experiment and signif- 
icantly simpler than the previous proposal fioj . The Ma- 
jorana fermions can be adiabatically moved around each 
other through control of the pinning potentials from fo- 
cused laser beams. We construct explicitly two composite 
braiding operations A and B, for which the different or- 
ders of braiding operations AB and BA, acting on the 
initial vacuum, give orthogonal output states that can 
be distinguished through the Raman spectroscopy. The 
proposed experiment thus provides a promising approach 
to directly probe the fundamental non-Abelian statistics 
of the anyons. Realization of Majorana fermions has 
been proposed in other systems fWL [ll| . however, for 
typical condensed matter materials, it is very challeng- 
ing to directly demonstrate the non-Abelian statistics of 
the Majorana fermions by braiding them in space. The 
proposed realization of Majorana fermions in an s-wave 
atomic superfiuid therefore offers a unique opportunity 
to adiabatically rotate them in space with focused laser 
beams to observe the non-Abelian fractional statistics, 
which is the most interesting property associated with 
these exotic quasi-particles. 

We consider a cloud of ultracold fermionic atoms with 
three relevant spin components in the ground-state man- 
ifold, label as |1),|2),|3) (e.g., they can be taken as 
the lowest three Zeeman levels of ^Li atoms near the 
broad s-wave Feshbach resonance). To generate the 
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effective SO coupling, we apply three traveling wave 
laser beams in the x and y directions which couple the 
atomic levels |1),|2),|3) resonantly to the excited state 
|e) through the standard tripod configuration as shown 
in Fig. 1 (l7l - [l9| . The single-particle Hamiltonian for 
each atom takes the form Hg — p^/2m + i?ej, where p 
denotes the momentum operator and m is the atomic 
mass. The light-atom interaction Hamiltonian Hf,j is 

given by H,, = ?iEj=?(f^j|e)(jl + h.c.) = hVi\e) {B\ + 
/i.e., where the corresponding Rabi frequencies take the 
form 0.1 — Vt sin/3 cos4>e~^^^, 0.2 = cos/3e~''^^, and 
5I3 = 51 sin/3 sin 06**''^ with k the wave number of the 
light, n = ^|fli|2 -t- |fl2p + If^ap, and the bright state 
\B) = (Sll\l) + Vll\2) + /n. We assume uniform 

plane wave laser beams so that fi,/3, are all constants 
and the values of /3, will be specified below. The di- 
agonalization of Hej yields two orthogonal dark states 
\Di) = sin(/.e-*'=(2'-^)|l) - cos(/ie-*'=(''+^)|3), ID2) = 
cos /3 cos (l)e-'^^y~'''> 1 1) - sin /3|2) + cos /3 sin ^e-*'=(2'+^) |3) , 
which are not coupled to the excited state |e) and thus 
immune to spontaneous emission. We assume the atoms 
are initially pumped to these dark states and they re- 
main in the dark states as long as the typical velocity v 
of the atoms satisfy the condition rjvk <C O. [i^l, where 
the prefactor 77 « 1.3 for the isotropic SO coupling dis- 
cussed below. In typical experiments, the wave number 
k 10^ m^^ and the Rabi frequency O of the order of 
magnitude 2tt x 20 MHz can be readily obtained, so the 
condition r]vk <^ V. is fully satisfied for ultracold atoms 
with the typical atomic velocity about several centime- 
ters per second. The probability of non-adiabatic excita- 
tion, estimated by {rjvk/V,)^, is negligible (about 10~^). 
In a resonant coupling configuration as shown in Fig. 
la, the effective atomic spontaneous emission rate is re- 
duced from the original rate 7^ by a factor of the non- 
adiabatic excitation probability and thus estimated by 
7s(?7wfc/ri)2, where is at a few Hz level. So the atomic 
spontaneous emissions induced by these laser beams are 
negligible if the experimental time is within 100 ms. The 
effective atomic spontaneous emission rate can be further 
reduced if one uses off-resonant coupling configuration or 
increases the Rabi frequency fl with more laser power. 
Under the above condition, we can project the Hamil- 
tonian Hs to the subspace spanned by these two dark 
states, and the projected Hamiltonian takes the form 
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where Ax = —hk cos (3ax, Ay = (?ifc/2) sin^ /3(Ty -t- 
{Tik/2) (1 -I- cos^ /3). In the derivation, we have assumed 
(j) — tt/4: and defined an effective spin with the basis 
vector I t) = + i\D2))/V2 and | i) = + 

\D2))/V2- The Pauli matrices axjCTy are defined in as- 
sociation with these spin basis vectors. Equation ^ de- 
scribes the atom in a spin-dependent synesthetic gauge 
field and a spin-independent version of this gauge 



field has been demonstrated in recent experiments 12 1| . 
To get the standard form of the SO coupling, we choose 
the value of /3 to satisfy sin^ (3 = 2 cos/? = 2{y/2 — 1) and 
observe the system in the moving frame with the veloc- 
ity Vy = {hk/2m) (1 -I- cos2/3), so that we can drop the 
constant term {hk/2) (1 -t- cos^ /3) in Ay. 

The effective Hamiltonian Hso has the required SO 
coupling, but we still need to generate a perpendicular 
Zeeman field on the effective spin | t) 1 1 i) with the cou- 
pling term hzaz. For this purpose, we apply two addi- 
tional laser beams that couple the states |1) and |3) to the 
excited state |e) off-resonantly with a large detuning 
(see Fig. 1(a)), with the corresponding Rabi frequencies 
denoted by fli^ = i |fl^| e"**^^ and = in^le''"'. With 
the driving \0.iz \ , \0.^z\ ^ ^d, it leads to a perturbation 
Hamiltonian given by — — i/jfipe^^'^^ll) (3| -I- /i.c, where 
Op = iJl^p/Ad |22|]. We assume Op < O, so the Hamil- 
tonian Hp can not pump the atoms outside of the dark 
state space spanned by the | f) and | |) levels, however, 
it splits the degeneracy of these two levels and leads to 
a coupling Hamiltonian of the form Hz ~ {hOp cos /3) az 
in the dark state space, which gives exactly the effective 
Zeeman field term hzaz with hz = hOp cos /3. The Zee- 
man field hz/h can be continuously tuned from zero to 
about 27T X 50 kHz (enough for our later applications), 
and under this condition, the requirement fip O is 
well satisfied with the off-resonant excitation probabil- 
ity {Op/0)'^ < 10"^. Note that in this scheme we need 
only two laser beams shined from the same direction as 
Oi,0.2 to create the required Zeeman field. This signifi- 
cantly simplifies the experimental requirement compared 
with the scheme in Ref. 19] which needs five additional 
laser beams and more complicated optical configuration 
to create the Zeeman field. 




FIG. 1; Schematic representation of the light-atom interac- 
tion configuration for generation of the effective Hamiltonian 
with SO coupling and a perpendicular Zeeman field, (a) The 
relevant atomic levels and the coupling Rabi frequencies to 
corresponding laser beams with specified polarizations a+, 
a-, or TT. (b) The spatial configuration and propagating di- 
rection of the laser beams. 

Now we have the effective single-particle Hamiltonian 
Hso + Hz of the right form. The atomic interaction 
is described by the s-wave scattering. In the second 
quantization formalism, it is described by the Hamil- 
tomanH.nt = E.< J W (r) (r) (r) d^r. 
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where (r) denote the fermionic field operators, i,j 
are summarized over the three spin components 1,2,3, 
and Qij are proportional to the s-wave scattering lengths 
in the corresponding channel. In the s-wave super- 
fluid state, the Hamiltonian Hint can be replaced by 
the elTective mean-field BCS type Hamiltonian Hbcs — 
J Ay (r) (r) (r) d^r + h.c, where (r) = 

9v correspond to the superfluid order pa- 

rameters. The typical s-wave pairing potential A^^ in 
experiments is about 2t: x 3.5 kHz [ij], which is much 
smaller than Q, so the interaction Hamiltonian Hbcs 
can not pump the atoms outside of the dark-state sub- 
space. Under this condition, we can project Hbcs to 
the dark-state subspace, and the projected Hamiltonian 
Hp takes the form Hp = J A (r) (r) vj/j^ (r) d^r -I- h.c, 
where A is a linear superposition of Ai2,Ai3,A23 
with the superposition coefficients determined from 
the transition matrix between |2), |3) and | f), 
I D, \B). If we write the wave- function ^(r) as 
$(r) — [u-f (r), u^(r), i;^(r), i;^(r)]^ in the basis speci- 
fied by the operators (^'t(r), 5';(r), ^'|(r), -5'|(r))^, the 
Schrodinger equation for the total Hamiltonian H = 
Hso + Hz + Hp with eigen-energy E takes the form of 
the Bogoliubov-de Gennes (BdG) equation 



Hso + Hz 



A(r) 



A*(r) -ay{Ht, + Ht)ay 



$(r) = E^{v). (2) 



If we introduce Tx^y^z to denote Pauli matrices in the 
Nambu spinor space, the Hamiltonian H can be written 
in the form 



H 



^i]Tz + a{cTxPx + 0-yPy)Tz + hzffz 



2™ 

-f (r^rR-eA - Tj,ImA) 



(3) 



where the laser induced SO coupling rate a — — ?ifc(^/2 — 
l). All the constant potentials have been included into 
the definition of the chemical potential fj, which is self- 
consistently determined from the particle density. For 
ultracold atoms in a global spin- independent confining 
trap V{r) , under the local density approximation, the 
global trap can also be included into the definition of the 
chemical potential with the substitution /i — !■ /i — V{r). 

The Hamiltonian H in Eq. ([3]) has been investigated in 
the context of a recent proposal for realization of non- 
Abelian anyons in a solid state heterostructure consist- 
ing of a superconductor, a 2D electron gas with Rashba 



SO coupling, and a magnetic insulator [11[. It is known 



that for this Hamiltonian, a Majorana fermion zero mode 
emerges from a vortex core if the parameters in the 

Hamiltonian satisfy the condition hz > he ~ y ^'^ + \ A\^ 



ll|,|l6|. For ultracold atoms, he is comparable with the 
Fermi energy Ef since the chemical potential \fj,\ ~ Ef 
and the pairing potential |A| ^ 0.35Ef. The condition 



hz > he can be easily satisfied since we assume hz/h can 
be tuned up to 2tt x 50 kHz, while Ef /his around 27r x 10 
kHz 1^. A vortex is needed for emergency of the Ma- 
jorana fermion in a central potential because the Hamil- 
tonian has the particle- hole symmetry ll|. Vortices in 
atomic superfluid can be created with several methods 
and have been observed in many experiments [23l 24 1. 



For instance, the required angular momentum associated 
with the vortices can come from initial rotation of the 
atomic cloud [23| or be imprinted from nonzero orbital 
angular momentum of laser beams 24 1 . 

We now turn to the important question of probing non- 
Abelian statistics associated with Majorana fermions. 
For an unambiguous demonstration of the non-Abelian 
statistics, we need to construct two braiding operations 
A and B. When applied to the initial state, the compos- 
ite braidings AB and BA give orthogonal output states 
which can be distinguished through the experimental de- 
tection. In the following, we explicitly construct these 
braiding operations for Majorana fermions and discuss 
the experimental detection of the output states. 

Let us consider 2A^ Majorana fermions confined in 
2N vortex cores, denote by 7i {i = 1, 2, • • • , 2iV). The 
Majorana fermions can be combined pairwise to form 
N complex fermionic states Cj = 72^ -I- 172^-1. Each 
fermionic state can be occupied or unoccupied and there- 
fore the ground state is 2^-fold degenerate. If we ex- 
change one Majorana fermion j with its nearest neigh- 
bor J -|- 1 by adiabatically moving the pinning poten- 
tial, we get an unitary operation Tj = exp (7r7j+i7j/4) — 
(l-|-7j+i7j)/\/2 acting on the 2^-dimensional degenerate 
Hilbert space [9] . The elementary interchange operations 
Tj {j — 1, 2, • • • , 2A^ — 1) form a complete set of genera- 
tors for the braid group. To demonstrate the non-Abelian 
statistics, we need at least four Majorana fermions due 
to existence of the superselection rule. For four Majo- 
rana fermions, they combine into two complex fermions 
ci,C2. The ground state has degeneracy four, but the 
Hamiltonian H conserves parity of the fermion number, 
and the even- number subspace spanned by {|0), cjcJlO)} 
is decoupled from the odd-number subspace spanned by 
{cj|0), C2IO)}. At low temperature (temperature signifi- 
cantly below the superfluid gap), the initial state is typ- 
ically a vacuum state |0) [12], and in the even- number 
subspace {|0), cjcjlO)}, only two of the three generators 
tit'''2, are independent with 

■ (4) 

From these generators, we construct two composite 
braiding operations 



A = TlT2T^ \ B = T^ ^T2T^ 



with the property 
AB 



IGz 



BA 



-IGx 



(5) 



(6) 
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where Tj~ represents a braiding along the reverse direc- 
tion. With the initial state |0), the braiding operations 
AB and BA yield two orthogonal output states, either 
i\0) or — ic|c2|0), demonstrating the non-Abelian nature 
of the anyons. During braiding, we need to fulfill the adi- 
abaticity. The gap to the first excited state in the vortex 
core is estimated by ujq ^ A'^/Ef ^ 2tt x 1.2 kHz. If the 
braiding time tb is of the order of tens of ms, the non- 
adiabatic excitation probability estimated by l/(aJoib)^ 
is negligible. 

The two output states i|0) and — icIcjiO) can be distin- 
guished through the Raman (or r.f.) spectroscopy [25| . 
For the state cjcJlO), we have two fermionic atoms at 
zero energy trapped at the vortex cores. We can apply 
a Raman 7r-pulse to resonantly bring these zero-energy 
atoms to a different hyperfine state \D) which is 

outside of Feshbach resonance and thus only weakly in- 
teracting with other atoms. If we choose the Rabi fre- 
quency of the Raman pulse to have a Gaussian shape 
= r2oexp(-r^/-u;^)exp(-Woi^) (where -tf < t < tf 
and the center of the vortex core is at r = 0), it will not 
break up the Cooper-pairs or excite other quasiparticles 
at the vortex core which are detuned typically by an en- 
ergy scale comparable with |A|. It is estimated that the 
error is about 4 x 10~^ and thus negligible under typical 
experimental parameters A = 27r x 3.5 kHz, ujq = 0.35A, 
ilo = 1-7^0, tf — S/wo and w — 1.5/zm. If the two laser 
beams for the Raman pulse propagate along different di- 
rections, it also applies a momentum kick to the target 
atoms to push them outside of the atomic cloud. These 
target atoms in a different hyperfine state \D) can then 
be detected with the standard quantum jump technique. 

In summary, we have proposed an experimental scheme 
to create and observe Majorana fermions in an s-wave 
superfluid of ultracold fermionic atoms with laser induced 
SO coupling. The non-Abelian statistics of the Majorana 
fermions can be directly probed through laser-controlled 
braiding and detection based on the Raman spectroscopy. 
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